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DESCRIPTION 

CARRIER PHASE GPS POSITIONING DEVICE AND METHOD 

5 TECHNICAL FIELD 

The present, invention relates to a carrier 
phase GPS (Global Positioning System) positioning device 
and method for detecting a position of a mobile station 
by carrier phase GPS positioning. 

10 

BACKGROUND ART 

Recently and continuing, carrier phase GPS 
positioning is widely used in the field of location 
survey . In the carrier phase GPS positioning, a receiver 

15 on a reference side and a receiver on a positioning side 
simultaneously receive signals from a plural number of 
satellites, and calculate accumulated values of the 
carrier phases of the satellite signals at the reference 
side and the positioning side, separately, resulting in a 

20 carrier phase accumulation value (abbreviated as "phase 
accumulation" if necessary, below) . The thus obtained 
carrier phase accumulation value contains an uncertainty 
factor that is integral multiple of the wavelength of the 
carrier. This uncertainty factor is referred to as 

25 "integral carrier phase ambiguity", and is often 
abbreviated as "integer ambiguity". 

A technique is well-known that uses a Kalman 
filter to determine the integer ambiguity. In this 
technique, a tracking filter is provided which regards 

30 the position to be determined and the integer ambiguity 
as state variables, a double phase difference of the 
phase accumulation on the positioning side relative to 
the reference side is an observation quantity, and each 
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time an observation is made the state variables are 
updated. 

There are also other techniques for 
determining the integer ambiguity. For example, it is 
5 known that the integer ambiguity related to a double 

phase difference can be found by the least-squares method 
under certain conditions by using the double phase 
difference of the carrier containing the integer 
ambiguity . 

10 In the related art, if the electromagnetic 

wave is interrupted after the integer ambiguity is 
determined or in the course of the determination (also 
referred to as * cycle-slip" ) , for example, the 
electromagnetic wave cannot be received, one has to 

15 determine the integer ambiguity again after reception of 
the electromagnetic wave resumes. However, because the 
above techniques of the related art are proposed 
specifically for positioning (determining the position 
of) an object which is fixed at a certain position for a 

20 long time, re-determination of the integer ambiguity is 
quite time-consuming . 

To solve this problem, there is a known 
technique in which, after reception of the 
electromagnetic wave resumes, a search space is 

25 established having a radius corresponding to a positional- 
variance with an output position of an IMU (Inertia 
Measurement Unit) as a center, and the integer ambiguity 
can be determined from the number of candidates of the 
integer ambiguity, which are solutions in the search 

30 space. For example, Japanese Laid on Patent Publication 
No. 2001-99919 discloses such a technique. 

However, in the above technique, in addition 
to a RTK positioning device, an IMU is also needed, and 
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one has to calculate the variance of the position 
measured by the IMU alone when the electromagnetic wave 
is interrupted. In addition, if the electromagnetic wave 
is cut off for a long time, the search space expends 
5 accordingly, and it is difficult to re-determine the 
integer ambiguity in a short duration. 



DISCLOSURE OF THE INVENTION 

Accordingly, it is a general object of the 

10 present invention to solve the above problem of the 
related art. 

A specific object of the present invention is 
to provide a carrier phase GPS positioning device capable 
of determining and re-determining an integer ambiguity 

15 quickly and precisely, a carrier phase GPS positioning 
method, a carrier phase GPS positioning system, and a. 
reference station . 

According to a first aspect of the present 
invention, there is provided a carrier phase GPS 

20 positioning device including a first integer ambiguity 

estimation unit that associates variance data in a first 
duration extracted from data received from a satellite by 
a reference station at a fixed position, with data 
received from the satellite by a mobile station in a 

2 5 second duration shorter than the first duration, and 

estimates an integer ambiguity included in a carrier 
phase accumulation value of signals transmitted from the 
satellite received by the mobile station; and a 
positioning unit that determines the position of the 

3 0 mobile station using the integer ambiguity estimated by 

the first integer ambiguity estimation unit. 

According to the present invention, because 
the integer ambiguity is estimated by using sampling data 
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on the mobile station side in the second duration shorter 
than the first duration, the time up to estimation of the 
integer ambiguity is shortened. The carrier phase GPS 
positioning device of the present invention can be 
5 implemented as a mobile station receiving data from a 
reference station, a reference station receiving data 
from a mobile station, or a device receiving data from 
both the reference station and the mobile station. 

It should be noted that the variance data and 

10 the data received by a mobile station are not necessarily 
to be of the same type. 

Preferably, abnormal values are excluded from 
the variance data. In addition, when reception of an 
electromagnetic wave emitted from the satellite is 

15 temporarily interrupted, data prior to the interruption 
is excluded from the variance data. 

Preferably, the variance data in the first 
duration includes a plurality of carrier phase 
accumulation values of the signals transmitted from the 

20 satellite at a first number of times in the first 

duration, and the first integer ambiguity estimation unit 
associates the plurality of carrier phase accumulation 
values on the reference station side at the first number 
of times, with a plurality of carrier phase accumulation 

25 values on the mobile station side at a second number of 

times, and estimates an integer ambiguity included in the 
carrier phase accumulation values of the signals 
transmitted from the satellite received by the mobile 
station, and here the second number is less than the 

30 first number. Alternatively, the first integer ambiguity 
estimation unit associates the plurality of carrier phase 
accumulation values on the reference station side at the 
first number of times, with a carrier phase accumulation 
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value on the mobile station side at one time, and 
estimates an integer ambiguity included in the carrier 
phase accumulation value of signals transmitted from the 
satellite received by the mobile station. In the latter 
5 case, single epoch positioning is possible. 

In addition, preferably, after the first 
integer ambiguity estimation unit estimates the integer 
ambiguity, the positioning unit determines the position 
of the mobile station using data measured on the mobile 

10 station side alone. 

Therefore, the amount of communication data 
between the mobile station and the reference station is 
greatly reduced after estimation of the integer ambiguity. 

In addition, preferably, the carrier phase GPS 

15 positioning device further includes a movement quantity 
detection unit that detects a movement of the mobile 
station and a movement quantity of the mobile station 
when the mobile station is moving, a second integer 
ambiguity estimation unit that, when the mobile station 

20 is at rest, estimates the integer ambiguity included in 
the carrier phase accumulation value of the signals 
transmitted from the satellite received by the mobile 
station. The estimation is made based on a plurality of 
carrier phase accumulation values on the reference 

25 station side and a plurality of carrier phase 

accumulation values on the mobile station side at a 
plurality of times in the period when the mobile station 
is at rest, and a third integer ambiguity estimation unit 
that, while the mobile station is moving, estimates the 

30 integer ambiguity included in the carrier phase 

accumulation value of the signals transmitted from the 
satellite received by the mobile station while taking 
movement detection results into consideration. 
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According to the present invention, the 
integer ambiguity estimation units carry out the 
estimation processing in parallel and independently from 
each other. Because integer ambiguities independent from 
5 each other are estimated , by comparing and investigating 
the integer ambiguities, an appropriate integer ambiguity 
can be obtained, and this can increase precision and 
reliability of the positioning. 

As an embodiment, after the second integer 

10 ambiguity estimation unit or the third integer ambiguity 
estimation unit estimates the integer ambiguity, the 
positioning unit determines the position of the mobile 
station using the integer ambiguity estimated by the 
second integer ambiguity estimation unit or the third 

15 integer ambiguity estimation unit instead of the integer 
ambiguity estimated by the first integer ambiguity 
estimation unit. 

Further, if the mobile station is a vehicle 
having wheels, the movement quantity detection unit 

20 detects a movement of the vehicle based on a wheel speed 
sensor that detects a rotational speed of the wheels. 
When a slip ratio greater than a predetermined value is 
detected by at least the wheel speed sensor, the integer 
ambiguity estimation processing by the third integer 

25 ambiguity estimation unit is initialized, and the 

positioning unit determines the position of the mobile 
station using the integer ambiguity estimated by the 
first integer ambiguity estimation unit until the third 
integer ambiguity estimation unit estimates or re- 

30 estimates the integer ambiguity. 

But when the integer ambiguity has been 
estimated by the second integer ambiguity estimation unit, 
the integer ambiguity estimated by the second integer 
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ambiguity estimation unit may be used for determining the 
position of the mobile station. 

Because the third integer ambiguity estimation 
unit takes movement quantity detection results into 
5 consideration, even when the mobile station is moving, it 
is possible to estimate the integer ambiguity at high 
precision . 

In addition, preferably, when plural reference 
stations are present in a communication region, a 

10 reference station is selected which is able to 

communicate with more satellites in common with the 
satellite communicating with the mobile station, and 
variance data related to the selected reference station 
are used. In addition, when there are plural reference 

15 stations able to communicate with the same number of 

satellites, a reference station is selected which has the 
highest minimum reception strength of signals from the 
satellites. Further, when plural reference stations, 
which receive signals from plural common satellites and 

20 the signal reception strength of each of the common 

satellites exceeds a predetermined value, are present in 
a communication region, a reference station is selected 
which is closest to the mobile station, and variance data 
related to the selected reference station are used. 

25 Therefore, even the reference station changes 

along with the movement of the mobile station, reduction 
of the integer ambiguity estimation precision is 
preventable . 

The carrier phase GPS positioning device of 

* 

30 the above inventions may be installed in the navigation 
device of a vehicle acting as a mobile station, or other 
movable objects such as a working robot, a mobile phone, 
and a PDA, or alternatively, in a facility capable of bi- 
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directional communication with the mobile station. 

According to a second aspect of the present 
invention, there is provided a carrier phase GPS 
positioning method, including the steps of associating 
5 variance data in a first duration extracted from data 

received by a reference station at a fixed position from 
a satellite, with data received from the satellite by a 
mobile station in a second duration shorter than the 
first duration, and estimating an integer ambiguity 

10 included in a carrier phase accumulation value of signals 
transmitted from the satellite received by the mobile 
station; and determining the position of the mobile 
station using the estimated integer ambiguity. 

According to a third aspect of the present 

15 invention, there is provided a carrier phase GPS 

positioning method including the steps of acquiring a 
carrier phase accumulation value at one time on the 
mobile station side; acquiring a plurality of carrier 
phase accumulation values at a plurality of times prior 

20 to the one time on the reference station side; 

associating the carrier phase accumulation values on the 
reference station side at the plural times, with a 
carrier phase accumulation value on the mobile station 
side at the one time, and estimating an integer ambiguity 

25 included in the carrier phase accumulation value of 

signals transmitted from the satellite received by the 
mobile station. 

According to a fourth aspect of the present 
invention, there is provided a carrier phase GPS 

30 positioning system including a reference station that 
extracts variance data in a first duration based on 
received data from a satellite; a carrier phase GPS 
positioning device including a first integer ambiguity 
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estimation unit that associates the variance data with 
data received from the satellite by a mobile station in a 
second duration shorter than the first duration, and 
estimates an integer ambiguity included in a carrier 
5 phase accumulation value of signals transmitted from the 
satellite received by the mobile station; a positioning 
unit that determines the position of the mobile station 
using the estimated integer ambiguity; and a 
communication path that enables communication between the 
10 carrier phase GPS positioning device and the reference 
station . 

According to a fifth aspect of the present 
invention, there is provided a reference station that 
extracts variance data in a predetermined duration based 

15 on received data from a satellite, and transmits the 

variance data to a carrier phase GPS positioning device 
including an estimation unit that associates the variance 
data with data received from the satellite by a mobile 
station, and estimates an integer ambiguity included in a 

20 carrier phase accumulation value received by the mobile 
station from the satellite; and a positioning unit that 
determines the position of the mobile station using the 
estimated integer ambiguity • 

According to a sixth aspect of the present 

25 invention, there is provided a reference station 

•including an acquisition unit that acquires a carrier 
phase accumulation value at one time on a mobile station 
side; an integer ambiguity estimation unit that 
associates a plurality of carrier phase accumulation 

30 values at a plurality of times prior to the one time on 
the reference station side, and estimates an integer 
ambiguity included in a carrier phase accumulation value 
of signals transmitted from the satellite received by the 
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mobile station; a positioning unit that determines the 
position of the mobile station using the integer 
ambiguity estimated by the integer ambiguity estimation 
unit; and a transmission unit that transmits the position 
5 detected by the positioning unit to the mobile station. 



BRIEF DESCRIPTION OF THE DRAWINGS 

These and other objects, features, and 
advantages of the present invention will become more 
10 apparent from the following detailed description of 
preferred embodiments given with reference to the 
accompanying drawings, in which: 

FIG. 1 is a schematic view of a carrier phase 
GPS positioning device according to the present 
1 5 invention ; 

FIG. 2 is a diagram showing a configuration of 
the carrier phase GPS positioning device in FIG. 1; 

FIG. 3 is a block diagram showing an 
embodiment of a carrier phase GPS positioning device 3 4 
2 0 installed in the mobile station 30 according to the 
present invention; 

FIG. 4 is a view illustrating the definitions 
of coordinate systems used in descriptions; 

FIG. 5 is a flowchart illustrating the method 
25 of determining the integer ambiguity in the carrier phase 
GPS positioning device 34 according to the present 
embodiment ; 

FIG. 6 is a flowchart illustrating the 
optional processing subsequent to the routine in FIG. 5; 
30 FIG. 7 is a flowchart illustrating an 

operation of determining the integer ambiguity by the 
time-series determination method of the related art, 
which is performed in parallel to the routine in FIG. 5 
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and/or the routine in FIG. 6; 

FIG. 8 is a flowchart illustrating an 
operation performed in the reference station 20 
corresponding to the operation performed in the mobile 
5 station 30 as shown in FIG. 5, FIG. 6, and FIG. 7; and 

FIG. 9 is a flowchart illustrating an 
operation of the mobile station 30 which is capable of 
communicating with a number of reference stations 20. 



10 BEST MODE FOR CARRYING OUT THE INVENTION 

Below, preferred embodiments of the present 
invention are explained with reference to the 
accompanying drawings . 

FIG. 1 is a schematic view of a carrier phase 
15 GPS positioning system according to the present invention. 

As illustrated in FIG. 1, the carrier phase 
GPS positioning system includes GPS satellites 10 
orbiting around the earth, a reference station 20 located 
at a fixed position (known position) , and a mobile 
20 station 30 that is on the earth, and is able to move on 
the earth. 

Each of the GPS satellites 10 broadcasts 
navigation messages to the earth continuously. The 
navigation messages include orbital information of the 

25 corresponding GPS satellite 10, clock correction value, 

and correction coefficients of the ionospheric layer. The 
navigation messages are spread using a C/A code, are 
carried on a LI carrier (frequency: 1575.42 MHz), and are 
broadcast toward the earth . 

30 Presently, there are 24 GPS satellites 

orbiting around the earth at an altitude of 20,000 km in 
the sky. Every four GPS satellites are equally arranged 
on one of six orbital planes of the earth, which are 
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inclined by 55 degrees relative to each other. Therefore, 
at least five satellites are always observable from a 
position as long as the position is open to the sky, no 
matter where the position is on the earth. 
5 FIG. 2 is a diagram showing a configuration of 

the carrier phase GPS positioning device in FIG. 1. 

In FIG. 2, the mobile station 30 has a GPS 
receiver 32. In the GPS receiver 32, there is an 
oscillator (not illustrated) having an oscillating 

10 frequency equal to the carrier frequency of the GPS 
satellite 10. The GPS receiver 32 converts an 
electromagnetic wave, which is emitted from the GPS 
satellite 10 and is received by the GPS receiver 32 via a 
GPS antenna 32a, and then performs C/A code 

15 synchronization using the C/A codes generated in the GPS 
receiver 32, and extracts the navigation messages. 

The GPS receiver 32 calculates a carrier phase 
accumulation value <E> iu of the carrier waves from the GPS 
satellites 10i. Here, in the phase accumulation value 3>i U , 

20 the subscript i (=1 , 2, „.) represents the numbers assigned 
to the GPS satellites 10i, and the subscript u represents 
that the accumulation value is calculated on the side of 
the mobile station 30. 

The phase accumulation value <E>± U can be 
25 described as the difference between a phase ®i U (t) of the 
oscillator at the time t of receiving the carrier wave 

and a phase ®i U (t-r) of the carrier wave when the 
satellite signal from the GPS satellite 10i is generated, 
as shown by the following formula (1) . 

30 

$iu(t)= ® iu (t) - ®iu(t-T u )+N iu +£ iu (t) (1) 



Here, t u represents travel time from the GPS 
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satellite 10 to the GPS receiver 32 , and e iu represents 
noise (uncertainty) . Further, at the time when starting 
observing the phase difference, the GPS receiver 32 can 
accurately determine the carrier phase within one 
5 wavelength of the carrier wave, but cannot determine what 
number of wavelengths the present wavelength is. For this 
reason, in the phase accumulation value ^i U (t), as shown 
in the formula (1) , there is an uncertainty factor N iu , 
known as "integer ambiguity". 

10 The mobile station 30 also includes a 

communication device 33, such as a mobile phone. As 
described below, the communication device 33 is capable 
of communicating with a communication facility 23 
installed on the reference station 20 side, such as a 

15 base station for mobile phones, by bi-directional 
communication . 

A GPS receiver 22 having a GPS antenna 22a is 
installed in the reference station 20. The GPS receiver 
22, the same as the GPS receiver 32 in the mobile station 

20 30, calculates a carrier phase accumulation value <3>ib at 
time t based on the carrier waves from the GPS satellites 
10±, as shown by the following formula (2) „ 



<E > ib(t)= ©ib(t)-© ib (t-T b )+N ib + e ^(t) (2) 

25 

Here, N ib is an integer ambiguity, and £ ib 
represents noise (uncertainty) . In the phase accumulation 

value ^ib, the subscripts b represents that the 
accumulation value is calculated on the side of the 
30 reference station 20. 

The reference station 20 transmits the 

obtained carrier phase accumulation value ^>±b to the 
mobile station 30 via the communication facility 23. More 
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than one reference stations 20 may be installed in a 
specified region. As illustrated in FIG. 2, each of the 
reference stations 20 may be connected to one or more 
communication facilities 23 through Internet or other 
5 networks, or a communication facility 23 may be installed 
in each of the reference stations 20. In the former case, 
as long as the mobile station 30 is able to communicate 
with the communication facility 23, the mobile station 30 
can obtain the information received by each of the 

10 reference stations 20. 

FIG. 3 is a block diagram showing an 
embodiment of a carrier phase GPS positioning device 3 4 
installed in the mobile station 30 according to the 
present invention . 

15 The carrier phase GPS positioning device 34 of 

the present embodiment includes a calculation unit 40, 
which is connected to the GPS receiver 32 and the 
communication device 33, and further, to various sensors 
50 in the mobile station 30. The calculation unit 40 may 

20 also be installed in the GPS receiver 32. When the mobile 
station is a vehicle, the GPS receiver 32, the 
calculation unit 40 and/or the communication device 33 
may also be mounted in a navigation device. 

The calculation unit 40 may be formed from a 

25 micro-computer, and may include a satellite position 
calculation unit 42 , a movement quantity introduction 
unit 44, a state variable extraction unit 46, and an 
integer ambiguity estimation unit 48, as illustrated in 
FIG. 3.- 

30 The satellite position calculation unit 42, 

based on the orbital information in the navigation 
messages received by the GPS receiver 32 , calculates 
positions (Xi(t), Yi (t) , Zi(t)) of all observable GPS 
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satellites 10i at time t in a world coordinate system. 

FIG. 4 is a view illustrating the definitions 
of coordinate systems used in the following description. 

FIG. 4 shows relations between the world 
5 coordinate system, a local coordinate system, and a body 
coordinate system . 

As illustrated in FIG. 4, in the world 
coordinate system, the origin is defined to be at the 
center of gravity of the earth, the X axis and Y axis, 

10 which are perpendicular to each other, are in the 

equatorial plane, and the Z axis is perpendicular to the 
X axis and Y axis. 

The body coordinate system is defined on the 
body of the vehicle.. 

15 Because movement of each of the GPS satellites 

10 is confined to an orbital plane passing through the 
center of gravity of the earth, and the orbit of each of 
the GPS satellites 10 is an ellipse with the center of 
gravity of the earth as a focus, positions of each of the 

20 GPS satellites 10 in the orbital plane can be calculated 
by successive numerical solutions of Kepler's equation. 

Because the orbital planes of each of the GPS 
satellites and the equatorial plane in the world 
coordinate system satisfy a rotational transformation 

25 relation, positions (Xi (t) , Yj. (t) , Zi (t) ) of the GPS 
satellites 10 at the time t of receiving the carrier 
waves can be calculated by three dimensional rotational 
coordinate transformation of the positions of the GPS 
satellites 10 on the orbital planes. 

30 The movement quantity introduction unit 44, 

based on the output signals of the various sensors 50 
input periodically, calculates quantities related to 
movement of the mobile station 30, and outputs the 
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results to the state variable extraction unit 46. 

For example, if the mobile station 30 is a 
vehicle, the movement quantity introduction unit 44 
calculates the speed Vx (t) (speed in the forward and 
5 backward direction) and Vy (t) (speed in the right and 
left direction) at the time t of receiving the carrier 
wave based on output signals from various sensors 50, for 
example, two wheel speed sensors mounted on the driven 
wheels of the vehicle, a yaw rate sensor, a left and 

10 right G acceleration sensor, an azimuth meter. 

Because the speed vector (Vx (t) , Vy (t) ) of 
the vehicle is defined in the body coordinate system, 
whose origin is on the body of the vehicle, it is 
necessary for the movement quantity introduction unit 44 

15 to transform the speed vector (Vx (t) , Vy (t) ) from the 

body coordinate system to the world coordinate system via 
the local coordinate system. Usually, the rotational 
transformation of coordinates can be performed by using 
Euler angles. In the present embodiment, the 

20 transformation from the body coordinate system to the 
local coordinate system is performed using only a yaw 

angle <l> (t) since the roll angle and pitch angle are 
small. Depending on the situation, the roll angle and 
pitch angle may also be considered, or the yaw angle may 
25 also be ignored. The transformation from the local 
coordinate system to the world coordinate system is 

performed by using the longitude <f> (t) and latitude 1 (t) 
of the position of the vehicle. 

Specifically, assume the position of the 
30 vehicle in the world coordinate system is (X u , Y u , Z u ) , 
and the longitude and latitude of the position of the 
vehicle is ( <t> , X ) , the speed vector of the vehicle in 
the world coordinate system d/dt[X u , Y u , Z u ] can be 



WO 2005/054889 



PCT/JP2004/018328 



- 17 - 



expressed by the following formula (3) 



d/dt[X u , Y u , Z u ] T =rot(<f> , X)*rot(<M*[Vx, Vy] 



(3) 



5 Here, [] T means transpose of a matrix, rot ( <t> , 

X ) and rot ( 0 ) are defined as the following formulae (4) 
and (5) . 



10 



— cos^sin^l 

— sin ^ sin /I 

cos A 



— sin <£ 
cosfi 

0 



(4) 



15. 



20 



25 



30 



rot(<p) - 



cos q> -smp 
sin <p cos <p 



(5) 



The longitude <|> (t) and latitude A, (t) of the 
position of the vehicle may be known (fixed longitude and 
latitude of a specified location which has been already 
determined, or may be separately measured variable 
longitude and latitude of the moving vehicle (that is, 
the mobile station 30) . 

The yaw angle <t> (t) may be calculated by 
integrating a yaw angular speed (the output signal of the 
yaw rate sensor) , or may be determined by using an 
azimuth meter. 

In the above formula (3) , by using input 
quantities U01 f U02, and U03 on the right side, and 
expressing the formula in a discrete way, one obtains the 
following formulae (6) , (7) , (8) . 



Xu(tn) - Xu(tn-i) + DT*U01 

Yu(tn) - Y u (tn-l) + DT*U02 
Zu(tn) - Zu(tn-l) + DT*U03 



(6) 
(7) 
(8) 
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Therefore, the final known input is expressed 

as below. 

5 U = [DT*U01, DT*U02, DT*U03] T (9) 

In formulae (6) , (7) , (8) , and (9) , DT 
represents a sample time interval (data updating 
interval) , and satisfies t n = t n -i + DT . Below , for 
10 simplicity of explanation, it is assumed that the sample 
time interval DT equals a period of calculating the phase 
accumulation value by the GPS receivers 22 and 32 . 

The integer ambiguity estimation unit 4 8 
estimates the integer ambiguity by using two types of 
15 carrier phase accumulation values, that is, the carrier 

phase accumulation value O ib on the side of the reference 
station 20, which is received by the mobile station 30 
through the communication device 33, and the carrier 

phase accumulation value $i u on the side of the mobile 
20 station 30. 

Specifically, a double phase difference of the 
phase accumulation of the GPS satellites 10j and 10h (j 
is not equal to h) at time t can be expressed by the 
following formula (10) . 

25 

*jhbu =(*jb(t)-* ju (t))- (*hb(t)-*hu(t>) (10) 

On the other hand, because the distance 
between the GPS satellite 10i and the GPS receiver 22 or 
30 32 equals the wavelength L of the carrier wave multiplied 
by the phase accumulation value, the double phase 

difference ^>jhbu of the phase accumulation satisfies the 
following formula (11) . 
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<J> .... = 



jhbu 




X b (.t)-Xj(t)y +(.Y b (t)-Y J (.t)y+(Z b (t)-Z J Q)f 
JiX u (r) - Xj (f)) 2 + (7 B (/) - Yj {t)f + (Z„ (* ) - Zj (* )) 2 } 



(11) 



(X b (t)-X h (t)) 2 +(Y b (t)-Y h (t)) 2 +(z,(*)-z A (0) 2 
V(X tt (0 - X, (0) 2 + (T. CO - r A (r)) 2 + (Z. W - Z h it)) 2 \ /L+N Jhhtt +e Jhbu 



In formula (11), [X b (t), Y b (t), Z b (t) ] are 
10 coordinates (known) of the reference station 2 0 at time t 
in the world coordinate system, and [X u (t), Y u (t), Z u (t) ] 
are coordinates (unknown) of the mobile station 30 at 
time t, tXj(t), Yj(t), Zj (t) ] and [X h (t), Y h (t), Z h (t)] are 
coordinates of the GPS satellites 10j and 10 h at time t 
15 calculated by the satellite position calculation unit 42. 
Njhbu represents the double phase difference of the 
integer ambiguity, that is, N jh bu = (N jb -N ju ) - (Nhb-N hu ) . 

The integer ambiguity estimation unit 4 8 

deduces double phase differences ^jhbu by using the 
20 formula (10) for each of four or more combinations of two 
GPS satellites 10j and 10 h , thereby obtaining relations 
satisfying the formula (11) by substituting the deduced 

double phase differences ^jhbu into the formula (11) . For 
example, when five GPS satellites 10i to 10s are 
25 observable, with the GPS satellite 10i as a reference 

satellite, ^12^, ^i3bu^ ^i4bu/ ^i5bu are substituted into 
the formula (11) and four relations are obtained. 

In this way, from the time of starting signal 
reception and periodically thereafter, that is, at time t 
30 = ti, t2, ... t n , the integer ambiguity estimation unit 4 8 

deduces the double phase differences ^jhbu by using the 
formula (10) and gives relations satisfying the formula 
(11) at each time of ti, t2, ... t n . Then, the integer 
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ambiguity estimation unit 48 calculates the variance of 

* jhbu on many data samples over a long time period (from 
ti to t n ) , for example, by the Kalman filter, least 
squares method or other estimation techniques , and 
5 determines integer ambiguities N ju , N hu (this processing 
is described below in detail with reference to FIG. 7) . 
In other words, the integer ambiguity estimation unit 48 

needs to store the double phase differences ^jhbu at each 
time before the integer ambiguities Nj U , N hu are 

10 determined. Further, the integer ambiguities Nj U , N hu 
related to the reference station 20 have already been 
determined on the side of the reference station 20, and 
are known quantities . The integer ambiguity estimation 
unit 48 acquires the integer ambiguities N ju , N hu related 

15 to the reference station 20 by communications. After the 
integer ambiguities Nj U , Nhu are determined, by well-known 
interf erometric positioning methods, accurate coordinates 
of the mobile station 30 can be deduced. 

In practical measurements, however, because of 

20 interruption of the electromagnetic wave, the 

communication between the GPS satellite 10i and the GPS 
receivers 22 or 32 may be interrupted sometimes. In this 
case, it is necessary to re-determine the integer 
ambiguity, and it is desired that the re-determination of 

25 the integer ambiguity be completed in a time period as 
short as possible. 

However, in the aforementioned method, many 
data samples over a long time period from ti to t n , are 
needed to determine the integer ambiguity, and hence the 

30 re-determination of the integer ambiguity is quite time- 
consuming, for example, 60 seconds or so. This problem 
occurs not only when communication with the GPS satellite 
is interrupted during measurements, but also at the time 
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when the measurement is started. 

In contrast, by the method of determining the 
integer ambiguity of the present embodiment, as described 
below, it is possible to re-determine the integer 
5 ambiguity in a short time. 

Below, the method of determining the integer 
ambiguity of the present embodiment is explained with 
reference to the accompanying drawings. For clarity, the 
aforementioned method is referred to as "time-series 

10 determination method" to be distinguished from the method 
of the present embodiment. 

FIG. 5 is a flowchart illustrating the method 
of determining the integer ambiguity in the carrier phase 
GPS positioning device 34 (specifically, the integer 

15 ambiguity estimation unit 48) according to the present 

embodiment. It is assumed that in addition to the routine 
in FIG. 5, the integer ambiguity estimation unit 48 is 
able to execute processing in the time-series 
determination method and processing of calculating the 

20 position of the mobile station 30 after the integer 
ambiguity is determined. 

In step S100, when the electromagnetic wave is 
interrupted, or when the ignition switch of the vehicle 
is switched ON, the routine in FIG. 5 is initiated. Here, 

25 "electromagnetic wave interruption' 7 means a state in 

which the GPS receiver 32 in the mobile station 30 can no 
longer receive signals required for positioning from the 
GPS satellites 10 (that is, the phase accumulation value 
becomes discontinuous, also referred to as "cycle-slip"). 

30 In this case, it is substantially not possible to execute 
processing in the aforesaid time-series determination 
method and processing of calculating the position of the 
mobile station 30 after the integer ambiguity is 
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determined . 

In step S105, if the electromagnetic wave is 
interrupted, the routine proceeds to step 110 after the 
reception of signals from the GPS satellites 10 is 
5 resumed (in this example , assume the timing of proceeding 
to step 110 is t = t n ) . 

If the ignition switch of the vehicle is 
switched ON, the routine proceeds to step 110 immediately. 

In step S110, the phase accumulations ^ibftn-a), 

10 ^ib(tn) prior to t = t n (here, from t n - a to t n ) , which 

are calculated on the side of the reference station 20 f 
are input to the integer ambiguity estimation unit 48. 
For example, these data may be acquired by sending an 
appropriate request signal from the mobile station 30. 

15 Alternatively, if the mobile station 30 is continuously 
receiving signals from the reference station 20, the 
integer ambiguity estimation unit 48 may read these data 
from a memory in the mobile station 30. In the latter 
case, the phase accumulation data prior to the present 

20 time, that is, 3>ib(tn-a), ^ibCtn-i) may also be read out. 

In the step S120, at the time t when the 
signal reception is resumed (in this example, assume t = 

t n ) , the phase accumulation value ^iu(t n ) on the side of 
the mobile station 30 is input to the integer ambiguity 

25 estimation unit 48. 

In the step S130, the integer ambiguity 
estimation unit 48 determines or re-determines the 
integer ambiguity based on the phase accumulation values 
obtained in the step S110 and S120. In this step, a 

3 0 method is employed to determine the integer ambiguity, 
which is different from the aforesaid time-series 
determination method in combinations of the phase 
accumulation values 9 used in deduction of the double 
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phase differences of the phase accumulation values. 

Specifically, in the step S130 f for example, 
the double phase differences O jhbu related to two GPS 
satellites 10j and 10 h (j is not equal to h) can be 
5 deduced by using the following formula (12) . 

^jhbu =(^jb(t k )-O ju (t n ) )- <*bb(t k )-* hu (t„) ) (12) 

In detail, in formula (12) , the phase 
10 accumulation values <I>ib (t n - a ) , O ib (t n ) , which are 

obtained in step S110, are substituted into $j b (t k ), 
respectively, and the phase accumulation values <I>hb(tn-a), 

^hb(tn) obtained in step S110 are substituted into <& 
jb (t k ), respectively. Thus, at time t=t n , a plural number 
15 of double phase differences <3> jh bu are generated (in this 
example, a number of a+1 double phase differences are 
generated). The *j U (t n ) and <E>hu(t n ) in formula (12) are 
the phase accumulation values$ ju (t n ) and <*>hu(t n ) at time 
t = t n obtained in step S120. 

2 0 Accordingly, a+1 equations are obtained by 

substitution into formula (11) . In formula (11) , the 
known coordinates of the reference station 20 at time t n - a 
to t n are substituted into [X b (t), Y b (t), Z b (t)], the 
known coordinates of the GPS satellites 10j and 10 h at 

2 5 time t n - a to t n are substituted into [Xj (t) , Yj (t) , Zj (t) ] 

and [X h (t) , Y h (t) , Z h (t)] with respect to the reference 
station 20, respectively, known coordinates of the GPS 
satellites 10j and 10 h at the time of t n (that is, [Xj (t n ) , 
Yj(tn), Zj(tn)] and [X h (t n ), Y h (t n ), Z h (t n )]) are always 

3 0 substituted into [Xj (t) , Yj (t) , Zj (t) ] and [X h (t) , Y h (t) , 

Zh(t)] with respect to the mobile station 30, 
respectively, and [X u (t), Y u (t), Z u (t)] are always the 
coordinates (unknown) of the mobile station 30 at t= t n . 
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These descriptions can be generally expressed 
by the following equation (13) . 



10 



jhbu 



= [{^x 6 a,)-x,a fc )) 2 +(y A (r it )-r y (r Jt )) 2 +(z A a A )-z y .(/,)) 2 
- > /(x 1( a n )-x y a„)) 2 +(r l( a n )-y,.(f B )) 2 +(^a n )-2 y a n )) 2 } 



(13) 



j lL+N jhbu +£ Jhbu 



Then, the integer ambiguity estimation unit 4 8 
calculates the variance of £ jhbu using these data, for 
example, by the Kalman filter, least squares method or 
15 other estimation techniques to determine the integer 
ambiguities Nj U , Nhu- 

It should be noted that in the method 
described in the present embodiment, the data after the 
time t=t n are not necessary, and the integer ambiguity 
20 can be determined at the time when reception of signals 

is resumed or started. Hence, single epoch positioning is 
possible. Due to this, it is possible to quickly 
(instantaneously) start or resume positioning at the time 
of starting the vehicle or after recovery from the cycle 
25 slip. 

In addition, because the integer ambiguity is 
an integer, in step S130, the integer ambiguity can be 
obtained by finding the integer solutions (that is, wave 
number) closest to the real solutions obtained by the 
30 Kalman filter. For example, the LA-MBDA can be used for 

this purpose, which un-correlates the integer ambiguities, 
and narrows the searching space of the integer solutions 
so as to facilitate finding the solutions. 
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After the integer ambiguities are determined, 
in. step S140, positioning is started or re-started by 
well-known interf erometric positioning methods (for 
example, kinematic positioning (RTK-GPS positioning 
5 algorithm) ) , and accurate coordinates of the mobile 
station 3 0 can be deduced. The position of the mobile 
station 3 0 obtained in this way may be used in various 
controls or be presented as information, for example, 
th.ey may be output and displayed on a screen of a 

10 navigation device, or displayed in a map shown in a 
screen of a mobile phone. 

As described above, according to the present 
embodiment , the integer ambiguity can be determined by 
combining data on the side of the mobile station 30 at a 

15 certain time and data on the side of the reference 

station 2 0 prior to the certain time. Therefore, it is 
possible to quickly (instantaneously) start or resume 
positioning at the time of starting the vehicle or after 
recovery from the cycle slip. 

2 0 In addition, in the routine illustrated in FIG. 

5, it is not always necessary to obtain all the phase 
accumulation values from t n - a to t n , from the reference 
station 20, but some of the phase accumulation values may 
be absent. In other words, as long as a sufficiently 

25 reliable integer ambiguity can be obtained in step S130, 
the number of the phase accumulation values obtained from 
the reference station 20 can be any number. 

In addition, in the routine illustrated in FIG. 
5, instead of obtaining the phase accumulation values 

30 from the reference station 20, variance data of the phase 
accumulation values may be obtained. For example, If the 
covariance matrix created by the integer ambiguity 
estimation unit 48 in step 130 can be created on the side 
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of the reference station 20 (for this purpose, it is 

necessary to transmit the phase accumulation value 9 
iu(t n ) to the reference station 20), the covariance matrix 
may be transmitted from the reference station 20 to the 
5 mobile station 30. From the same point of view, the 

mobile station 30 may transmit the covariance matrix to 
-the reference station 20, the processing in step 130 and 
position calculation of the mobile station 30 may be 
executed on the side of the reference station 20, and the 

10 reference station 20 sends only the integer ambiguity 
and/or the position of the mobile station 30 to the 
mobile station 30. In this case, the processing load of 
"the mobile station 30 can be greatly reduced, and the 
amount of communication data between the reference 

15 station 20 and the mobile station 30 can be greatly 
reduced. 

In addition, in the routine illustrated in FIG. 
5, it is exemplified that there is only one phase 

accumulation value (^iuftn)) on the side of the mobile - 
20 station 30, which is used for determining the integer 
ambiguity, but the number of the phase accumulation 
values on the side of the mobile station 30 may be more 
"than one, and two or more phase accumulation values (less 
■than a+1 in the above example) can be used in the same 
25 manner. For example, when using data after the time t = 
-tn, although the time required for determining the 
integer ambiguity increases, many combinations can be 
obtained even with a small number of the phase 
accumulation values obtained from the reference station 
30 20, and similarly, it is possible to determine a 
sufficiently reliable integer ambiguity. 

In addition, in the routine illustrated in FIG. 
5, for example, when five GPS satellites 10i to 10 5 are 
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observable, and the cycle-slip occurs in the GPS 
satellite 10 2 , it is not necessary to carry out the 
processing of re-determination of the integer ambiguity 
for GPS satellites other than GPS satellite 10 2 . In this 
5 case, for example, for the double phase difference <I>i2bu 
of the phase accumulation values, the aforementioned 
formulae (12) and (13) can be used by using <3>iu(t n ), ® 
2u(t n ). For the other double phase differences ($ 13bu , <I> 
i4bu, ^lsbu) of the phase accumulation values, the 
10 aforementioned formulae (10) and (11) can be used by 

using <£ lu (t n - a ) ^iu(tn) r 0 3u (t n - a ) , ... ^ 3 u(t n ). In this 
case, estimation of the integer ambiguities N ju , N hu is 
performed by using a combination of these data. 

FIG. 6 is a flowchart illustrating the 

15 optional processing subsequent to step S140 in FIG. 5 
(that is, after the time t = t n ) . For example, the 
routine in FIG. 6 may be executed when the phase 
accumulation value <J> iu (t n ) obtained in step S120, or the 
position of the mobile station 30 deduced based on the 

20 phase accumulation value ^iu(t n ) are apparently 

unreasonable, for example, the change of the position of 
the mobile station 30 during the cycle-slip is not 
practical according to the moving speed of the mobile 
station 30 . 

25 In step S150, at the sample time t=t n +i next to 

the time t = t n , the phase accumulation value ^i U (t n +i), 
which is calculated on the side of the mobile station 30, 
is input to the integer ambiguity estimation unit 48. It 
should be noted that after the time t = tn, if the phase 

30 accumulation value on the side of the mobile station 30 
cannot be observed again because of the cycle slip, the 
routine in FIG. 5 is carried out again. 

In step SI 60, the phase accumulations O ib (t n -a- 
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1) , »w ^ibCtn+i) prior to t = t n +i (here, from t n - a -i to t n +i) , 
which are calculated on the side of the reference station 
20, are input to the integer ambiguity estimation unit 48. 
For example, these data may be acquired by sending an 
5 appropriate request signal from the mobile station 30. 
Alternatively, if the mobile station 30 is continuously 
receiving signals from the reference station 20, the 
integer ambiguity estimation unit 48 may read these data 
from a memory in the mobile station 30. 
10 In step S170, the same processing as that in 

the step S130 in FIG. 5 is performed by using the phase 

accumulation value ^iu(t n +i) . Namely, in formulae (12) and 
(13), t n is replaced with t n +i • The integer ambiguity 
estimation unit 48 processes these data and determines 

15 the integer ambiguities N ju , N hu . In addition, in step 

S170, estimation of the integer ambiguities N ju , N hu may 
be performed by combining the phase accumulation values 
obtained in the present step and the phase accumulation 
values obtained in the previous steps. For example, if 

20 the present routine is executed for the first time, the 
phase accumulation values obtained previously are the 
phase accumulation values obtained in step S130 in FIG. 5. 
If the present routine is not executed for the first time, 
the phase accumulation values obtained previously are the 

25 phase accumulation values obtained in step S130 in FIG. 5 
and the phase accumulation values obtained in the 
previous step S170. 

In step S180, the integer ambiguity previously 
determined is compared with the integer ambiguity 

30 presently obtained. For example, if it is determined that 
the integer: ambiguity obtained in step S130 is obviously 
abnormal fxrom history of the speed of the mobile station 
30 during the cycle-slip, positioning can be continued 
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with the integer ambiguity presently obtained (refer to 
step S140) . 

In step S190, step S150 to step S180 are 
repeatedly executed until an integer ambiguity is 
5 determined by the aforesaid time-series determination 

method, which is executed at the same time as the present 
routine (time subsequent to t=t n , that is, t=t n +m, (m=2 , 3 , 
...) ) . 

If an integer ambiguity is determined by the 
10 time-series determination method, in step S200, the 

integer ambiguities determined with the two methods are 
compared with each other, and then positioning is 
continued with the integer ambiguity determined by the 
time-series determination method (that is, taking over of 
15 the integer ambiguity can be performed) by well-known 
interf erometric positioning methods. 

FIG. 7 is a flowchart illustrating an 
operation of determining the integer ambiguity by the 
time-series determination method, which is performed in 
2 0 parallel to the routine in FIG. 5 and/or the routine in 
FIG. 6. 

In step S300, it is determined whether the 
mobile station 30 is moving or not. This determination 
may be executed by the wheel speed sensor, when the 
25 mobile station 30 is a vehicle, or image processing. 

* 

If it is determined that the mobile station 30 
is at rest, the routine proceeds to step S310. 

If it is determined that the mobile station 30 
is moving, the routine proceeds to step S340. 
30 In step S310, a static model is used to 

execute the time-series determination method to determine 
the .integer ambiguity. In step S310, if the mobile 
station 30 starts to move, the processing is stopped, and 
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initialized. On "the other hand, if the mobile station 30 
is at rest until the integer ambiguity is determined f the 
routine proceeds to step S370. In step S370, the value of 
the thus determined integer ambiguity is compared with 
5 the value of the integer ambiguity obtained by the 
routines in FIG. 5 and FIG. 6 (refer to step S200) . 

In step S340, a movement model is used to 
execute the time— series determination method to determine 
the integer ambiguity. Specifically, the integer 
10 ambiguity estimation unit 48 establishes the following 

state equation indicated by (14) by using the known input 
introduced by the movement quantity introduction unit 44 
(refer to formula. (9)) . 

15 V (t„)= V (tn-l) + U(tn-l) + W(ta-l) (14) 

Here, v (t n ) is a state variable at t = tn, and 
may represent position coordinates (unknown) of the 
mobile station 30 [X u (t n ) , Y u (t n ) , Z u (t n )], and the double 
20 phase difference of the integer ambiguity Njhbu- U and W 
are the aforesaid known input and external noise (system 
noise: normal white noise), respectively. In the static 
model, the item of the known input U(t n _i) in formula (14) 
does not exist. 

25 In addition, the integer ambiguity estimation 

unit 48 also establishes the following observation 
equation indicated by (15) , which holds also in the 
static model. 

30 Z(t n )= H(t n )*77(t n ) + V(t„) (15) 

Here, Z and V represent an observation 
quantity and observation noise (normal white noise) , 
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respectively. The observation quantity Z is the double 
phase difference of the integer ambiguity N jh bu (refer to 
formula (10)) . While the state equation (14) is a linear 
equation , the observation quantity Z is non-linear 
5 relative to the state variable X u/ Y u , Z u , the items in 

the formula (11) are partially differentiated relative to 
Xu / Y u , Z u , thereby obtaining H in formula (15) . 

Therefore, if applying the Kalman filter to 
the state equation (14) and the observation equation (15) , 
10 the following equations can be obtained. 

For updating the time, 

V (t n ) (_> = 7} (t n -i) ( + > + U (t n -l) + W (t n -l) (16) 
P(t n ) M - P(t n -i) (+) + Q(tn-l) (17) 



15 



For updating the observation, 



K(t n )=P(t n _i) ( " } *H T (t a ) * (H(t n ) *P (t n ) M *H T ( t„ ) +R ( t„) ) ~* (18) 

V (t n ) <+) = V (t n ) (_) + K(tn)*(Z(tn)-H(t n )*77 (t n ) M ) (19) 
20 P(t n ) (+> =P(t n ) ( - ) - K(t n ) *H(t n ) *P (t n ) ( " > (20) 



Here, Q and R represent the covariance matrix 
of the external noise and the covariance matrix of the 
observation noise, respectively. The formulae (13) and 

25 (19) are covariance equations. Here, the superscript ( " } 
and <+) indicate time before and after the updating, 
respectively. In the static model, the item U(t n -i) in 
formula (16) does not exist. 

As a. result, the estimated values of the 

30 integer ambiguity are found to be real solutions. However, 
since the integer ambiguity is actually an integer, the 
integer ambiguity is found to be an integer solution 
(that is, wave number) closest to the real solutions. For 
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example, the LA-MBDA can be used for this purpose, which 
un-correlates the integer ambiguities , narrows the 
searching space of the integer solutions so as to 
facilitate finding the solutions . 
5 In step S350, if it is determined the slip 

ratio of the mobile station 30 exceeds a preset value, 
since the reliability of the movement model (that is, the 
aforesaid known input) declines when the slip ratio of 
the mobile station 30 exceeds a preset value, the 

10 processing is stopped, and initialized in step S360. When 
the mobile station 30 is a vehicle, the slip ratio of the 
mobile station 30 can be detected based on an output 
signal from the wheel speed sensor or accelerator sensor, 
as is well known in the field of the ABS control. 

15 If it is determined that the slip ratio of the 

mobile station 30 is does not exceed a preset value, the 
routine proceeds to step 370. 

In step 370, the value of the thus determined 
integer ambiguity is compared with the value of the 

20 integer ambiguity obtained by the routines in FIG. 5 and ■ 
FIG. 6 (refer to step S200) . 

After the processing is stopped as shown in 
step S330 and S360, the integer ambiguity determined by 
the routines in FIG. 5 or FIG. 6 is used until the 

25 integer ambiguity is determined by either of the models. 

In the processing shown in FIG. 7, when cycle 
slip occurs, the processing is stopped and initialized, 
similar to steps S330 and S360, and afterward, the 
integer ambiguity determined by the routines in FIG. 5 or 

30 FIG. 6 is used until the integer ambiguity is determined 
by either of the models. 

As described above, according to the present 
embodiment, by introducing a movement model (that is, 
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input a known external input U(t) to the Kalman filter), 
the integer ambiguity can be determined precisely even 
when the mobile station 30 is moving. In addition, by 
separately using the movement model and the static model 
5 depending on the movement condition of the mobile station 
30, it is possible to improve the estimation precision of 
the integer ambiguity. By setting the known input to the 
movement model to be zerro when the mobile station 30 is 
at rest, the movement model alone is also sufficient. 

10 In addition, £>y executing the processing in 

FIG. 7 in parallel to the routine in FIG. 5 and the 
routine in FIG. 6, the two independent integer 
ambiguities can be compared with each other, and an 
appropriate integer ambiguity can be selected to be used 

15 in positioning, hence enhancing positioning reliability. 
From this point of view, when the GPS receivers 22 and 32 
are two-frequency receivers capable of receiving both of 
the LI wave and L2 wave emitted from the GPS satellite 10, 
for each of the LI wave and L2 wave, the same estimation 

20 processing can be performed in parallel simultaneously. 
In this case, because th.e sum of the periods of the two 
sides (Wide-Lane) can be created, and thereby, it is 
possible to further narrow the range of the candidates of 
the integer solutions. 

25 FIG. 8 is a flowchart illustrating an 

operation performed in the reference station 2 0 
corresponding to the operation performed in the mobile 
station 30 as shown in FIG. 5, FIG. 6, and FIG. 7. 

In step S400 , as a routine process, each of a 

30 number of reference stations 20 stores and manages the 
phase accumulation values O ib , which are measured based 
on the satellite signals from each of the GPS satellites 
10i, in correspondence to each of sample times and each 
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of the GPS satellites 10±. This treatment may be executed 
by a center facility (not illustrated) which controls the 
GPS satellites 10i. 

In step 410, as long as the cycle slip does 
5 not occur concerning data related to a GPS satellite 10k, 
the processing in step 400 is performed continuously. Of 
course, each time new data are created, data prior to a 
certain time can be deleted sequentially. In addition, 
integer ambiguity Nib between each of the reference 
10 stations 20 and each of the GPS satellites 10i can be 
introduced into processing of step 400, which is at an 
initial stage. Each of the reference stations 2 0 may 
determine the reliability of the obtained phase 

accumulation value <I>ib, and, for example, if abnormal 

15 values exist, it may be required not to store the 

abnormal values, or not to transmit the abnormal values 
to the mobile station 30. 

The data stored and extracted in step S400 may 
be transmitted to the mobile station 30 via the 

20 communication facility 23 upon receiving a request signal 
from the mobile station 30, as described in step S110. 
Alternatively, the data may be broadcast in each cycle of 
receiving signals from each of the GPS satellites 10±. 

In the former case, determination of a mobile 

25 station 30 as a transmission destination is accomplished 
by an ID code included in the request signal from the 
mobile station 30. In the latter case, the reference 
station 20 may send only special data to the mobile 
station 30 in response to the request from the mobile 

30 station 30. 

In step S400, if it is detected that the cycle 
slip occurs concerning a GPS satellite 10k, the routine 
proceeds to step S4 20, and the reference station 20 
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erases data stored before the cycle slip related to the 
GPS satellite 10k occurs. 

As described above, according to the present 
embodiment, it is not necessary to transmit data 
5 generated before the cycle slip occurs to the mobile 
station 30. Alternatively, the stored data generated 
before the cycle slip occurs may also be held, but not to 
be transmitted to the mobile station 30, that is, only 
data generated after the cycle slip occurs are 

10 transmitted to the mobile station 30. Thus, the mobile 
station 30 can obtain highly reliable data from the 
reference station 20 by just transmitting a request 
signal, and does not need to check the reliability of the 
data transmitted from the reference station 20. 

15 FIG. 9 is a flow chart illustrating an 

operation of the mobile station 30 which is capable of 
communicating with a number of reference stations 20. 

In step S500, it is determined whether the 
mobile station 30 is able to communicate with more than 

20 one reference stations 20. 

If it is determined that the mobile station 30 
can communicate with more than, one reference stations 20, 
the routine proceeds to step S510. 

In step S510, the correspondence relation 

25 between the GPS satellite 10i captured by the mobile 
station 30 and the GPS satellite 10± captured by one 
reference station 20 is detected. At the same time, the 
number of the GPS satellites 10, from which both the 
mobile station 30 and the reference stations 20 receive 

30 signals (below referred to as "common satellite 10") , is 
determined. 

In step S520, if it is determined that there 
is one reference station 20 which receives signals from 
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the largest number of common satellites 10, the mobile 
station 3 0 adopts the data from the one reference station 
20 and executes the routines illustrated in FIG. 5, FIG . 
6, and FIG. 7. For example , considering the case in which 
5 the mobile station 30 can receive signals from seven GPS 
satellites 10i to IO7, one reference station 20 x receives 
signals from six GPS satellites 10i to 10 6 , and another 
one reference station 20 2 receives signals from GPS 
satellites 10 x to 10 5 , and GPS satellites 10 8 to 10n, 

10 because both the mobile station 30 and reference stations 
20 receive signals from GPS satellites 10i to 10 6 , GPS 
satellites 10i to 10 6 are the common satellites 10, and 
because in this example, ttie reference station 20i 
receives signals from the largest number (six) of common 

15 GPS satellites 10i to 10 6 , the reference station 20i is 
selected. 

On the other hand, if it is determined that 
there are more than one reference stations 20 each of 
which receives signals front* the largest number of common 

20 satellites 10, the routine proceeds to step 530. 

In step 530, the reception strength from each 
of the common satellites 10 at different reference 
stations 20 are compared with each other. The reception 
strength at different reference stations 20 can be 

25 measured at the reference stations 20, and may be 
transmitted to the mobile station 30. 

In step S540, if it is determined that there 
is only one reference station 20 which receives signals 
from the common satellites 10 at a reception level higher 

30 than a reference level, the mobile station 30 adopts the 
data from the one reference station 20 and executes the 
routines illustrated in FIG. 5, FIG. 6, FIG. 7. In this 
case, a reference station may be selected which has the 
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highest minimum reception strength of signals from the 
common satellites. In this case, the subsequent step 550 
can be ignored. 

On the other hand, if it is determined that 
there are more than one reference stations 2 0 each of 
which receive signals from the common satellites 10 at 
the best reception condition, the routine proceeds to 
step S550. 

In step 550, distances to the mobile station 
30 from different reference stations 20 are compared with 
each other. 

In step 560, the reference station 20 having 
the shortest distance to the mobile station 30 is 
selected. As a result, the mobile station 30 adopts the 
15 data from the closest reference station 20, and executes 
the routines illustrated in FIG. 5, FIG. 6, and FIG. 7. 

As described above, according to the present 
embodiment, because data from the most appropriate 
reference station 20 are used among- the common satellites 
20 10 with which the mobile station 30 is able to 

communicate, it is possible to prevent degradation of the 
integer ambiguity estimation precision and the position 
detection precision when the reference station 20 is 
changed along with the movement of the mobile station 30. 
25 The priority order of selecting the reference 

station 20 is not limited to the above examples (namely, 
the number of the common satellites 10, reception 
strength, distance) , for example, under the condition 
that the number of the common satellites 10 is larger 
than a specified number and the reception strength 
exceeds a minimum level, the reference station 20 closest 
to the mobile station 30 may be preferentially selected. 
In this case (that is, distance is preferential), it is 



30 
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possible. to reduce the influence of the ionospheric layer 
refraction effect and tropospheric bending, and improve 
the estimation accuracy of the integer ambiguity. 

While the present invention has been described 
5 with reference to specific embodiments chosen for purpose 
of illustration, it should be apparent that the invention 
is not limited to these embodiments, but numerous 
modifications could be made thereto by those skilled in 
the art without departing from the basxc concept and 

10 scope of the invention. 

For example, in the above embodiments, it is 
described that the Kalman filter is applied to the state 
equation (14) and the observation equation (15) , but the 
least-squares or other estimation methods can also be 

15 used for estimating state quantities. 

In addition, in the above embodiments, it is 
described that the double phase difference is calculated 
to eliminate influence of the initial phase of 
oscillators in the GPS receivers 22 and. 32, and clock 

20 uncertainties. Of course, single phase difference can 
also be used as long as the influence of the initial 
phase of oscillators in the GPS receivers 22 and 32, and 
clock uncertainties can be eliminated. In addition, in 
the above embodiments, the influence of the ionospheric 

25 layer refraction effect, tropospheric bending effect, and 
multi-path are not considered, but the present invention 
is also applicable when considering these effects. 

In the above, for simplicity r the GPS 
satellite 10i is regarded as a reference satellite, but 

30 depending on the positions of the mobile station 30 and 

the reference station 20, any other GPS satellites 10 may 
also be used as the reference satellite. In addition, as 
long as there are four or more double phase differences 
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which are related to GPS satellites common to the mobile 
station 30 and the reference station 20 , any combination 
of the GPS satellites can be used for calculating the 
double phase differences. 
5 In the above embodiment, a vehicle is raised 

as an example of the mobile station 30. The mobile 
station 30 may also include a folk lift or a robot with 
the receiver 32 and/or the calculation unit 40, and a 
mobile 

10 mobile phone or PDA with the receiver 32 and/or the 
calculation unit 40. 

According to the carrier phase GPS positioning 
device of the present invention, it is possible to 
determine an integer ambiguity quickly and precisely. 

15 This patent application is based on Japanese 

Priority Patent Application No. 200 3-403640 filed on 
December, 2, 2003, the entire contents of which are 
hereby incorporated by reference. 



